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Abstract 

We study the rates allowed for the Higgs-mediated decays B^^ ^ fir, er and r — > 
/i/z/i, e/z/i in supersymmetric seesaw models, assuming that the only source of lepton 
flavour violation (LFV) is the renormalization of soft supersymmetry-breaking terms 
due to off-diagonal singlet-neutrino Yukawa interactions. These decays are strongly 
correlated with, and constrained by, the branching ratios for ^ — > /x/x and r 
/i(e)7. Parametrizing the singlet-neutrino Yukawa couplings and masses M^r- in 
terms of low-energy neutrino data, and allowing the flavour-universal soft masses 
for sleptons and for squarks, as well as those for the two Higgs doublets, to be 
different at the unification scale, we scan systematically over the model parameter 
space. Neutrino data and the present experimental constraints set upper limits on 
the Higgs-mediated LFV decay rates Br{B^ fJ-Tjer) ^4 x 10~^ and Br{T — > 
/i/x/x, e/x/i) < 4 X 10"^°. 
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1. Introduction. In the minimal supersymmetric extension of the Standard Model 
(MSSM), non-holomorphic Yukawa interactions of the form D^QH^ are generated at one- 
loop level [Q. At large values of tan/5 = V2/V1, the contribution of such loop-suppressed 
operators to the down-type quark masses may become comparable to those from the usual 
superpotential terms V^QHi. As these two contributions have different flavour structures, 
because of the up-type quark Yukawa interactions, they cannot be diagonalized simultane- 
ously 0, H . This difference leads to potentially large new contributions to Higgs-mediated 
flavour- changing processes involving down-type quarks, such as 5°^ — /i/i @, |^, 0, §] 
and mixing H, I . 

This discussion cannot be generalized directly to the charged-lepton sector, because 
the MSSM has no right-handed neutrinos at the electroweak scale. However, experimental 
data convincingly indicate that neutrinos do have masses 0. Their small values are most 
naturally explained via the seesaw mechanism [^, which involves super-heavy singlet 
(right-handed) neutrinos with masses M^r^. In this case, the presence of neutrino Yukawa 
couplings N^{Yy)ijLjH2 above the heavy-neutrino decoupling scale induces off-diagonal 
elements in the left-slepton mass matrix via renormalization [|ll|, 0, 



even if the initial slepton soft supersymmetry-breaking masses mo and trilinear couplings 
Aq are flavour- universal at Mqut- This is the only source of lepton-flavour violation (LFV) 
in supersymmetric seesaw models with flavour-universal soft mass terms. Since it is induced 
by heavy-neutrino Yukawa interactions, it relates the LFV processes to low-energy neutrino 
data. At the one-loop level, (|l|) also gives rise to flavour violation in non-holomorphic 
interactions of the form E'^LH2 and leads to Higgs-mediated LFV processes in the charged- 
lepton sector. 

It is well known that, at large tan/3, new Higgs-mediated contributions to the decay 
— > fifi may exceed the Standard Model (SM) contribution by orders of magnitude. 
The dominant contributions come from the diagrams in Fig. |I] (a) which are determined 
by Cabibbo-Kobayashi-Maskawa (CKM) matrix elements if the squark mass matrices are 
flavour- universal at the GUT scale, as suggested by the data on flavour- changing neutral 
interactions [|13| . It was suggested recently |]T3[ that the Higgs-mediated contribution to the 
LFV process r fififi depicted in Fig. |l| (b) might be sizeable in supersymmetric seesaw 
models Q, with a branching ratio as large as Br{T — >■ /x/x/x) ~ (9(10~^). This claim was 
made without a complete study of all related LFV processes. 

Since the only source of LFV in this model is (0), and all LFV processes are induced 
by slepton-neutralino and sneutrino-chargino loops, Higgs-mediated LFV is constrained 
indirectly by limits on the processes r /X7, r — 67 and fi 67, which have been well 
studied by experiment. Although r /i(e)7 and Higgs-mediated LFV depend on the model 
parameters in systematically different ways, the former impose important and unavoidable 
constraints on the latter processes. Additional constraints are expected to come from 
lepton-flavour-conserving but quark-flavour-violating processes such as —>■ fifi, which 

^ For a discussion of LFV decays in non-supersymmetric models the reader is referred to ^ . 
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Figure 1: Dominant Higgs penguin diagrams contributing to (a) ^ —>■ fifi and (b) t 
fififi decays at large tan/3. 
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Figure 2: Dominant Higgs penguin diagrams contributing to B^^ 
tan p. 
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depends on Higgs boson masses in the same way as Higgs- mediated r — > fifxfi (see Fig. jT]). 
In addition, there are other Higgs-mediated LFV processes such as B^ /ir and 5° — > er, 
which are induced by double-penguin diagrams, as depicted in Fig. 0, and may also be of 
experimental interest. The current experimental upper limits on such decays are shown 
in Table |l]. Although the double-penguin diagram is suppressed by an additional loop 
factor, its amplitude is enhanced by rrir/mf^ tan /5, as seen by comparing Fig. g with Fig. 0. 
Therefore, at large tan /3, the rate for 5° /ir may even exceed that of r — > yU/i/x. 

To quantify these qualitative statements, we adopt for the moment the approximation 

0, 



considered in 



z.e., we take all the supersymmetry-breaking mass parameters of the 
model to be equal at low scales, use heavy-neutrino masses that are degenerate with M^r = 
10^^ GeV, and assume that (1^^1^1^)32,33 = 1- This approximation is not realistic, but it is 
useful for comparing the sensitivities of different processes to new physics: a more complete 
treatment of these processes is presented in the following Sections of this Letter. In the 
simplified case, working in the mass- insertion approximation and following exactly we 
obtain 



Br{T 3fi) ~ 1.6 X 10" 



tan f5 
~60 



1 6 



100 GeV 



1 4 



(2) 



which is about six times smaller than the estimate quoted in that paper. This should be 
compared with the corresponding estimate 



Br{r n-f) ~ 1.3 x 10 



-3 
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60 



2 r 



100 GeV 



H 4 
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Both equations (0) and (^ are valid in the large-tan /3 limit. Whereas @ is two orders of 
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Table 1: Current experimental bounds for the branching ratios of the leptonic B decays 
Bs d — > ith ■ 



magnitude below the present experimental bound on r ^ 3/i, (Q) is three orders of magni- 
tude above the present bound on r — >• /i7. There is also the photonic penguin contribution 
to the decay r — > 3/i, which is related to Brij fij) by 



Br{T - 

Numerically, (||) leads to 
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Br{T 3/i)^ ~ 3.0 X 10^ 



tan (3 
60 



100 GeV 



(5) 



which is a factor of 100 larger than (|^). Notice also that suppressing @ by postulating 
large slepton masses would suppress (0) at the same time, since sleptons enter into both 
loops. However, suppressing (||) by a large Higgs mass Ma would not affect r — /i7. 

With the same assumptions, we obtain 



Br{Bl /i/i) ~ 1.9 X 10" 



tan/5 

"eo" 



n 6 



100 GeV 



Ma 



where only the leading tan j] dependence is presented, and 



Br{B^^ Tfi) ~ 3.6 X 10"^ 



tan P 
60 



100 GeV 

Ma 



(6) 



(7) 



to be compared with the upper limits in Table |^. In the case of Bd mesons, one should 
just multiply with |Vtd/Vtsp ~ 0.05. As expected, the Higgs-mediated branching ratio 
for fi° — > TyU can be larger than the one for r — > 3^. 

We comment in passing that the decays Bd^s /^e are suppressed by the ratio m?^/m1 ~ 
0.0036 compared to Bd^s ~^ f^T and, moreover, they are strongly constrained by the process 
/i — i> 67. We find that their branching ratios are below the range of prospective experimental 
interest. 
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2. Effective Lagrangians and branching ratios for LFV processes. We now present the 
calculational details we use in arriving at the approximate results of the previous Section 
and the numerical results to be presented in the next Section. We consider the i?-parity 
conserving superpotential: 



(8) 



where the indices i,j run over three generations and Mn is the heavy singlet-neutrino mass 
matrix. We work in a basis where (Ydjij, (Ye)ij and {MN)ij are real and diagonal. At the 



(9) 



one-loop level, this leads to the effective Lagrangian 0, 

-C^^f = 4Frf45,,ff? + (eo5i,+ey(FjFj,,)i/°*]4 + /i.c.+ 



where eo, ey, ei, and €2 are loop-induced form factors, and Eij = {Am^)ij/ml is the unique 
source of LFV. In the following, we present analytical expressions in the mass-insertion 
approximation, which is simple and known to reproduce well the full results in the case of 
large tan j3. However, in our numerical analyses of the next Section we use the full diagram- 
matic calculation of the (6s)-Higgs transition of [Q, which we modify according to to 
include the resummed large-tan (3 contributions to the down-type Yukawa couplings. Also, 
the LFV matrix Eij is calculated exactly using the numerical solutions to the full set of 



renormalization-group equations of the MSSM, including singlet neutrinos [|T2 
In the mass-insertion approximation, the down-quark form factors are 
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In the lepton sector, the corresponding form factors are 
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The resulting effective Lagrangians describing ffavour- violating neutral Higgs interactions 



with down quarks and charged leptons are [g, |I4 



reff 



eff 



(2Gl)^/^^^^^^ (diR d,L) [cos(a - + sin(a - - + h.c. , (15) 
cos"' p 

(2G'^)i/4!^ (r.i? /,x) [cos(a -I3)h^ + sin(a - /J)//" - lA^] + /i.e. , (16) 

COS kJ 



where 



[1 + (eo + eyYt^Su) tan/3] [1 + eotan/3] 



[1 + (ei + es^ii) tan/5] 
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(17) 
(18) 



As we are interested in S-meson decays, we have A^^ = V^lVtq, where the Vij are CKM 
matrix elements. 

Using (p!6[), one easily obtains the branching ratio for the decay r fififi, given by 
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(19) 



where Tr is the r lifetime and the large-tan /5 limit is taken in the last step. This is the 
result that was used to derive the estimate (§). 

The branching ratio for r — /i7 in the mass-insertion approximation reads ]T^ : 

a 



Br{T /i7) = ^mlTr\AL\'^ , 



(20) 



where 
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where D[f{x); x]{xi, X2) = (/(xi) -/(x2))/(xi -X2), and 
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This is the resuh that was used to derive the estimate (||). 

The dominant operators at large tan jS in the effective Hamiltonian describing the tran- 
sition b —>■ sltl~ are 



H 



GlMl 



w 



+ h.C. 



(22) 



where V is the CKM matrix and Gp the Fermi couphng constant. Here cs^p^io and Os,p,io 
are, respectively, the Wilson coefficients and the local operators 0, which are given by 



mbihRSL)ladj ■ 



(23) 



In the SM with the seesaw mechanism, the only non-zero coefficients are the cio- For « 7^ j, 
they scale like the square of the inverse mass of the singlet neutrinos, and are completely 
negligible. In supersymmetric seesaw models, there are two additional operators in (p2D, the 
scalar Os and the pseudoscalar Op, whose coefficients dominate over the SM contributions. 
In particular, for i ^ j they are suppressed by the scale of supersymmetry-breaking soft 
masses only. At large tan (3 the dominant contributions to and Cp come from the penguin 
and double-penguin diagrams presented in Fig. |I] and Fig. ||, respectively. We should note 
here, however, that when the Higgs masses are (9(TeV), and the sneutrinos are very light, 
the box diagrams may dominate. In this case we obtain in general small branching ratios, 
Br\Bs fir) ^ 10~^^, when we take the t ^ fi'j constraint into account. 

By setting cio to zero in (^), we obtain the following branching ratio 0, 



Br(fi, 
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JBs^bA ■ I 
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(24) 



where Mp^ and Tp^ are the mass and lifetime of the Bg meson, and fp^ = 230 ± 30 GeV 
is its decay constant. 

For Bg — ^ /i/x decay, the form factors are P, ||: 



d 



GpM^ cos^pX' 



bs 



sin(a — P) cos a cos(a — P) sin a 



Ml, 



Ml. 



^In (22) we have omitted operators proportional to the strange-quark mass, since they are subleading 
for our processes. 

•^We note that can be straightforwardly extended to Bd LFV decays by replacing s d. 
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and, for the double-penguin contribution to Bg fir decay, we obtain from ([TsD and ([T6|) 
the form factors 



T/l 



GfM^ cos^(3X 



sin^(Q; — /5) _|_ cos^(q; 
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The last two results are those 
(0). In the calculation for the 
the operators {hPL{R)s){jiPL{R)T). 



that were 
ratio Bs - 
However, 



used to derive the 
fiT one also has 
their contribution 



estimates and 
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which vanishes approximately at large tan /3. Furthermore, the 

operator {bPRs){flPLT) is proportional to m^m^ and thus subdominant to {hPLs){jlPjiT) 
we consider here. For the same reason, the leptonic CP-asymmetries in B^ fir decays 
due to the complex in (^) are of order m^/rriT-. 

3. Numerical Analyses. The purpose of this work is to study in a complete way the 
allowed rates for the Higgs-mediated LFV processes in supersymmetric seesaw models in 
which the only source of LFV is the renormalization of the soft supersymmetry-breaking 
mass parameters above Mtv^, due to the singlet-neutrino Yukawa couplings. We work 
with two models. First., we study the constrained MSSM (CMSSM) which has just two 
universal mass parameters at GUT scale, M1/2 for gauginos and mo for all scalars, taking 
for simplicity Aq = Q and not expecting our results to be sensitive to its exact value. 
Secondly., we study the more general flavour-universal MSSM (GFU-MSSM) in which the 
universal masses for squarks, sleptons and the Higgs doublets Hi and H2 are different 
from each other. This permits different mass scales for squarks and sleptons which, in 
turn, are independent of the Higgs boson masses. However, we always require that squark 
and slepton mass matrices at the GUT scale are each proportional to unit matrices. If 
one goes beyond this assumption, arbitrary sources of flavour violation appear in the soft 
supersymmetry-breaking sector, and the model loses all the predictivity, in particular the 
connection between the LFV and the neutrino masses and mixings. Moreover, LFV rates 

m- 

in terms of 



generally exceed experimental limits 



We parametrize the singlet-neutrino Yukawa couplings Yy and masses Mp^. 
low energy neutrino data according to We generate all the free parameters of the model 
randomly and calculate the low-energy sparticle masses and mixings by solving numerically 
the one-loop renormalizat ion-group equations and imposing the requirement of electroweak 
symmetry breaking. We calculate the rates for the LFV processes as described in the last 
Section. For the decays /j Ijj, we use the exact diagrammatic formulae in 



We fix the light-neutrino parameters by [^] Amig = 3 x 10"'^ eV^, Am^i = 5 x 10~^ 
eV^, tan^ 6'23 = 1, tan^ 6^12 = 0.4 and sin6'i3 = 0.1, corresponding to the LMA solution for 
the solar neutrino anomaly. The neutrino mixing phase is taken to be maximal: 6 = 7t/2. 
The overall light-neutrino mass scale is randomly generated in the range (10~^ — * 10~^) 
eV, with a constant logarithmic measure. We assume the normal mass hierarchy, which is 



favoured by leptogenesis arguments p3| . The rest of the neutrino parameters are calculated 



from the randomly generated textures Hi (which maximizes r — /i violation) and H2 (which 



maximizes r — e violation), that were introduced in ||21 



We fix tan (3 = 60, which is the largest value for which electroweak symmetry breaking 
is obtained generically. In the case of the CMSSM, the mass parameters M1/2 and mo for 
sleptons, squarks and both Higgses are generated randomly in the range (0 — > 700) GeV. 
In the case of the GFU-MSSM, we have four scalar mass parameters m^, mg, m^^ and 
for squarks, sleptons and Higgs doublets, respectively. We generate each of them randomly 
and independently in the range (0 700) GeV. We impose the experimental upper bounds 
Br{T ^ /i7) < 2 X 10"^ Br{T ^ 67) < 2 x 10"^ and 5r(/i 67) < 1.2 x 10~^\ and 
the lower bounds on unobserved particle masses are taken to be mj, rriq, Ma > 100 GeV, 
and no other phenomenological bounds are imposed. Thus, our analysis should give an 
accurate comparison of different LFV decays and Bg^d — ^ /W/x with each other, but may be 
somewhat optimistic concerning the overall rates. 

In Fig. we present scatter plots of Higgs-mediated i?r(r 3/i) against i?r(r fi'j) 
in (a) the CMSSM and (b) the GFU-MSSM. Whilst 5r(r ^ S^u) < 10"^^ in the CMSSM, 
in the GFU-MSSM the present 5r(r fi'j) bound allows 5r(r ^ 3/i) < 4 x 10~^°. This 
increase could have been expected, since in the GFU-MSSM the Higgs masses and slepton 
masses are independent of each other. The points with the highest possible Br{T — > 3/i) 
in In Fig. ^ correspond to Ma ~ 100 GeV. However, for tan/5 = 60, such low values 
for Ma are constrained by Bg /i/i. To see this, we plot in Fig. § Br^r 3/i) against 
Br{Bs —>■ fifi) in both models we consider. It follows that, in the GFU-MSSM, the present 
bound Br{Bs —>■ /i/i) < 2 x 10~^ constrains i?r(r —>■ 3/i) < 1 x 10~^°. This should be 
compared with the rate for the photonic penguin contribution to i?r(r 3/i), which is 
given by (|^). Our results show that, after imposing the t ^ fi'j and Bg fifi bounds, the 
photonic penguin dominates over the Higgs-mediated contribution. 

We emphasize that we have been seeking to maximize LFV effects in these figures. We 
also note that the lower bounds on the LFV processes in the plots are artificial, being due 
to the lower bounds on the generated values of Yi,. 

Proceeding with the S- meson decays, we plot in Fig. ^the Higgs-mediated Br{Bs /ir) 
against Br{Bs fifi) in (a) the CMSSM and (b) the GFU-MSSM. Again, there is a sharp 
upper bound on Br{Bs fir) in both models. In the GFU-MSSM the maximum allowed 
value for Br{Bs ^t) is 4 x 10~^, which is larger than that for r 3/i. However, detecting 
r leptons is experimentally challenging, particularly in a high-rate environment such as the 
LHC. To see the dependence of Bs — > /xr on Ma, we plot in Fig. |^ the value of Br{Bs — > /ir) 
against Ma- We emphasize that Ma can be very small in the GFU-MSSM |^^, whilst in 
the CMSSM there is an approximate lower bound Ma > 200 GeV, because it is related 
to the slepton masses. The points with the largest Br{Bs —>■ fir) in Fig. pi as well as the 
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points with the largest Br{r 3fi) in Fig. ||, correspond to Ma = O(IOO) GeV. 

Similar results are also valid for the decays r — > efifi and Bg —>■ er. Since phenomeno- 



logically the r — e transition can be as large as the t — one we have the same 



bounds Br{T —>■ efifi) < 4 x 10^^'^ and Br{Bs —>■ er) < 4 x 10~^. The scatter plots for 
those processes are indistinguishable from Figs. |^ to so we do not present them here. We 
have also studied the decays /le and /i eee. The Higgs-mediated contributions 

to Br{B^g^ — > /ie) and i?r(/i — >• eee) are suppressed by the — >• e7 constraint and small 
Yukawa couplings to be below 10~^^ and 10~^^, respectively. 

4- Conclusions. 

In this Letter we have studied the allowed rates for Higgs-mediated LFV decays in super- 
symmetric seesaw models, assuming that LFV is generated entirely by the renormalization 
effects of neutrino Yukawa couplings. Even if we allow the mass scales for squarks, sleptons 
and two Higgs doublets to differ from each other, the bounds due to the decays r fij, 
T —>■ ej, fi ^ e'-f and Bg —>■ fifj^ are very constraining. We obtain the following bounds on the 
Higgs-mediates processes: Br{Bg — > fir, er) ^ 4 x 10~^ and i?r(r — > /i/i/x, e/x/i) ^ 4 x 10~^°. 
We have discussed these numerical results in terms of approximate analytical expressions, 
but our calculations are more exact. We conclude that the Higgs-mediated contributions 
to r — > /i/i/i, efifi are subleading compared to the photonic penguin ones. 
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Figure 3: Scatter plot of Higgs-mediated Br{T — > S/i) against Br{T — > /I'j) in (a) 
CMSSM and (b) the GFU-MSSM. 
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Figure 4: Scatter plot of Higgs-mediated Br{T — >• 3//.) against Br{Bs — >■ /x/^) m ('a^ 
CMSSM and (b) the GFU-MSSM. 
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Figure 5: Scatter plot of Higgs-mediated Br{Bs — > //r) against Br{Bs — > m (a) the 
CMSSM and (h) the GFU-MSSM. 
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Figure 6: Scatter plot of Higgs-mediated Br{Bs /it) against the pseudoscalar Higgs mass 
Ma m (a) the CMSSM and (h) the GFU-MSSM. 
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